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Abstract 


Due to environmental and fuel cost concerns more and more wind- and solar- 
based distributed generation (DG) units are embedded in distribution networks 
(DNs). It is, however, a well-known fact that renewable energy generators are 
highly fluctuating sources, and therefore, energy storage systems such as battery 
storage systems (BSSs) are considered as a solution to handle such fluctuations. 
In general, DG units and/or BSSs convert traditional passive DNs (PDNs) into 
active DNs (ADNs). Consequently, it is important to investigate the impact and 
benefits of integrating such entities in conventional DNs. 

This dissertation presents a systematic study consisting of modeling, 
simulation, and optimization of dynamic operations of energy supply networks 
with embedded renewable generation and storage. Based on complex power 
flow models, different optimization problems are mathematically formulated 
and solved. 

In this work, novel mathematical models and a new combined problem 
formulation for active-reactive optimal power flow (A-R-OPF) in PDNs 
(without DG units and BSSs) and ADNs (with DG units and BSSs) are studied. 
Typically, DNs consist of two different networks in terms of voltage levels, 
namely, low-voltage and medium-voltage DNs. For this reason, investigations 
are carried out separately on both networks. Modeling procedures for PDNs, 
ADNs, and energy prices are presented, These procedures serve as the basis for 
this work. Then, simulation studies in PDNs are made to analyze its operating 
characteristics, In particular, the operation of on-load-tap-changers of main 
transformers is highlighted. Moreover, an optimization framework is introduced 


to minimize the total energy losses in PDNs. 


In ADNs, two voltage levels with two real case studies are separately 
considered. On the low-voltage level, a high penetration level of photovoltaic 
(PV) systems (PVSs) is considered in the network in order to reveal the impact 
of such a scenario. In particular, the reactive power capability of the inverters of 
these PVSs is explored. The total revenue from the installed PVSs is maximized 
whilst the total cost of energy losses and demand is minimized. Using different 
price models many interesting results are found, e.g., no need to use BSSs in 
low-voltage DNs for accommodating expected spilled PV energy. On the 
medium-voltage level, a DN with a high penetration of wind energy and BSSs is 
considered. In this case, the total revenue from wind parks and BSSs is 
maximized and the total cost of energy losses is minimized. It is found that a 
huge reduction in energy losses and reactive energy imports can be achieved. To 
prolong the life of BSSs only one fixed charge/discharge cycle every day is 
considered. The solution provides an optimal operation strategy which ensures 
the feasibility and enhances the revenue significantly. However, due to the fact 
that the profiles of renewable energy generation, demand and prices vary from 
day to day a fixed operation of BSSs cannot be optimal. 

A flexible battery management system is proposed to adapt to such 
variations. This is accomplished by optimizing the lengths (hours) of charge and 
discharge periods of BSSs for each day, leading to a complex mixed-integer 
nonlinear program (MINLP), An iterative two-stage framework is proposed to 
address this problem, In the upper stage, the integer variables (i.e., hours of 
charge and discharge periods) are optimized and delivered to the lower stage. In 
the Jower stage the A-R-OPF problem is solved by a NLP solver and the 
resulting objective function value is brought to the upper stage for the next 
iteration, This procedure will converge when number of iterations is reached. 


Using this flexible system a considerably higher revenue can be achieved. 
IV 


Zusammenfassung 


Bedingt durch Umweltbelange und steigende Kosten fiir fossile Brennstoffe 
werden immer mehr Wind- und Solaranlagen (distributed generation, DG) in 
Verteilernetzen (distribution networks, DNs) installiert. Es ist eine bekannte 
Tatsache, dass die Einspeisung durch erneuerbare Energietriger starken 
Schwankungen unterliegt. Ein méglicher Lésungsansatz zur Behandlung dieser 
Schwankungen ist die Nutzung von Energiespeichersystemen wie z. B. 
Batteriespeichersysteme (BSS). Der Einsatz solcher Systeme verwandelt 
traditionelle passive Verteilernetze (PDNs) in aktive Verteilernetze (ADNs). 
Folglich ist es wichtig, die Auswirkungen und Vorteile der Integration solcher 
Einheiten in konventionelle Verteilernetze zu untersuchen. 

In dieser Dissertation wird eine systematische Untersuchung (bestehend aus 
Modellierung, Simulation und Optimierung) des dynamischen Betriebs von 
Energieversorgungsnetzen mit eingebetteten erneuerbaren Energietragern und 
Speichersystemen vorgenommen. Basierend auf komplexen Lastflussmodellen 
werden verschiedene Optimierungsprobleme mathematisch formuliert und 
gelést. 

In dieser Arbeit werden neue mathematische Modelle und eine neue 
Problemformulierung fiir den kombinierten optimalen Lastfluss von Wirk- und 
Blindleistung (active-reactive optimal power flow, A-R-OPF) in PDNs (ohne 
DG-Anlagen und BSSs) und ADNs (mit DG-Anlagen und BSSs) vorgestellt. 
Typischerweise enthallen DNs zwei Spannungsebenen, nimlich Nieder- und 
Mittelspannung. Deshalb werden Untersuchungen in den _ beiden 
Spannungsebenen getrennt durchgeftihrt. Modellierungsverfahren fiir PDNs, 
ADNs und Energiepreise werden vorgestellt. Diese Verfahren dienen als 


Grundlage der vorliegenden Arbeit, Darauf  aufbauend — werden 
V 


Simulationsstudien in’ PDNs zur Analyse der  Betricbscigenschaften 
durehgetithrt. Insbesondere wird die Auswirkung des Betrichs der 
Laststufenschalter der Haupttransformatoren hervorgehoben. Dariiber hinaus 
wird ein Optimierungsverfahren zur Minimierung der gesamten Energieverluste 
in PDNs vorgestellt. 

In ADNs werden zwei Spannungsebenen mit jeweils zugehérigen realen 
Fallstudien getrennt betrachtet. Auf der Niederspannungsebene wird eine hohe 
Einspeisungsrate von PV-Anlagen (photovoltaic systems, PVSs) angenommen, 
um die Auswirkungen eines solchen Szenarios zu zeigen. Insbesondere wird die 
Fahigkeit der Inverter dieser PV-Anlagen zur Erzeugung von Blindleistung 
untersucht. Die Gesamteinnahmen aus den installierten PV-Anlagen werden 
maximiert, wahrend gleichzeitig die Gesamtkosten der Energieverluste und die 
Nachfrage minimiert werden. Durch die Verwendung unterschiedlicher 
Preismodelle kénnen viele interessante Ergebnisse generiert werden, z. B. 
besteht keine Notwendigkeit, BSSs in Niederspannungsnetzen fiir die Aufnahme 
iiberschiissiger PV-Energie zu installieren. Auf der Mittelspannungsebene wird 
ein DN mit einer hohen Einspeisungsrate von Windenergie und BSSs betrachtet. 
In diesem Fall werden die Gesamteinnahmen der Windparks und BSSs 
maximiert, wahrend die Gesamtkosten der Energieverluste minimiert werden. Es 
zeigt sich, dass eine enorme Reduktion der Energieverluste und der 
Blindleistungsimporte erreicht werden kann. Um die Lebensdauer der BSSs_ zu 
verlangern wird nur ein fester Lade-/ Entlade-Zyklus pro Tag betrachtet. Diese 
Lésung liefert eine optimale Betriebsstrategie, welche die Zuliissigkeit 
gewdhrleistet und den Profit signifikant erhdht. Aufgrund der Tatsache, dass die 
Profile der erneuerbaren Energien, der Nachfrage und der Preise von Tag zu Tag 


variieren, ist cin feststehender Betrieb der BSSs allerdings nicht optimal, 


Vi 


Weiterhin wird ein flexibles Batterie-Management-System zur Behandlung 
solcher Schwankungen vorgestellt. Dies wird durch die Optimierung der Lade- 
und Entladezeiten der BSSs fiir jeden Tag erreicht. Daraus resultiert ein 
komplexes gemischt-ganzzahliges nichtlineares Optimierungsproblem (mixed- 
integer nonlinear program, MINLP). Fiir dessen Lésung wird ein iteratives 
zweistufiges Verfahren eingefiihrt. In der oberen Stufe werden die ganzzahligen 
Variablen (d. h. Lade- und Entladezeiten) optimiert und an die untere Stufe 
weitergegeben. In der unteren Stufe wird das A-R-OPF Problem mit einem 
NLP-Léser gelést und der resultierende Wert der Zielfunktion wird an die obere 
Stufe fiir die nachste Iteration weitergegeben. Dieses Verfahren konvergiert, 
wenn eine Anzahl von Iterationen erreicht ist. Die Verwendung dieses flexiblen 


Ansatzes resultiert in bedeutend héheren Profiten. 
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Appendix A 


Appendix A: IEEE-RTS load data 


Table A.1: Load data of the low- and medium-voltage DNs (Hourly demand as a percentage 
of the annual peak demand) [50][7] 
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Appendix B 


Appendix B: Data for the low-voltage DN 


Table B.1: Data of the low-voltage DN [4] 


No. From Length R, Xx 
ae ae us [Ba us (km) (ohm/km) (ohm/km) 


ieee a 
eee 0157 [100 | e100 
ees ics | 1.900 | 0.100 


5 
Pp 5 of 2 | 6 | 0107 | 1900 | 0100 
ee ee 02 | 1900 0.100 
|. Bi 
eee 0.081 | 0383 | 
Po [| 9 | tw | =60070 6 S868 | 0078 
| io | io | i | 0093 | 086s | 0.078 
eee 010 | 0.195 0.070 
eemrmmee 0073 | i740 
ee eee 0119 | sig | 0350 
meen. 014s |i? | 04105 | 
eae 0035 | 900 | 0.100 
p28] 8 29 | 0084900 | 0.100 
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Table B.2: Data of the low-voltage DN (Demand daily peak and PFs) [41] 


Bus Preasli)(KW) Bus Pyeali)(KW) 


Table B.3: Data of the low-voltage ADN (Upper bounds of active and reactive power in 
forward and reverse direction at slack bus) [41] 


Table B.4: Data of the low-voltage ADN (PVSs) [41] 


at Specs max.pv(i) (kVA) = Ppy(i) Spcs.maxpv(i) (kVA) i Ppy(i) 
(kW) (kW) 
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Appendix C: Data for the medium-voltage DN 


Table C.1: Data of the medium-voltage DN [5] 


From | To Length R; xX B; 
ie ae Bus (km) (ohm/km) | (ohm/km) (us/km 


p i | 2 | 65.7000 | 0.169111 | 0.418206 | 3.9540 _| 
ee epee) 10100 | osteo -catszoe | 3.0500 — 
pS | 5 | 6 | 0.1300 | Ot69iit | 0.418206 | 3.9540 _| 
| 6 | 5 | 7 {| 0.1700 | Ot69111 | 0.418206 | 3.9540 _| 
| 7 | 7 | 9 | 0.2600 | 0.169111 | 0.418206 | 3.9540 _| 
| 8 | 9 {| 10 | 0.1400 | O169111 | 0.418206 | 3.9540 _| 
| 9 | 9 | ii | 0.3800 | Ot69111 | 0.418206 | 3.9540 _| 
p10 | it | 12 | 05600 | O.169111 | 0.418206 | 3.9540 _| 
| 14 | io | 17 | 1.0800 | 0.169111 | 0.418206 | 3.9540 _| 
| 16 | is | 19 | 0.4700 | 0.169111 | 0.418206 | 3.9540 _| 
| 17 | 19 | 20 | 0.4700 | 0.348124 | 0.468482 | 3.7571 _ 
p19 | 19 | 23 | 0.1900 | 0.348124 | 0.468482 | 3.7571 
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Table C.2: Data of the medium-voltage DN (Demand daily peak and PFs) [40] 


0.641346 0.95 


| 0.641346 | 0.95 | 25 | 0.028975 | 0.95 
02042 
.010¢ 


F 
0.001900 0.020425 


0.010450 
0.216600 


Table C.3: Data of the medium-voltage ADN (Upper bounds of active and reactive power in 
forward and reverse direction at slack bus) [40] 


Ssimax = 20 MVA 


Table C.4: Data of the medium-voltage ADN (wind turbines, PCSs capabilities and BSSs 
capacities) [40] 


Table C.5: Data of the medium-voltage DN (active energy prices for 24-hour-tariff price 
model in winter and spring) [42] 


Price Price 

h ($/MWh) ($/MWh) po ($/MWh) 
| 55.65 | 4643 | 9 | 78.91 | 70.02 | 17 | 82.23 | 66.33 _| 
| 52.33 | 45.70 | 10 | 79.47 | 72.97 | 18 | 83.07 | 67.81 | 
| 49.84 | 44.22 [it | 79.47 | 73.71 | 19 | 83.07_| 70.76 _| 


7 
71.44 | 62.65 | 16 | 78.08 24 52.33 51.59 
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Appendix D: Software implementation of DSI-1 


Here is the implementation of the DSI-1 used for carrying out dynamic 
power flow studies in DNs. The DSI-1 is implemented in the MATLAB- 
Simulink environment with user-interfaces, as shown in Figs. D.1 and D.2. 


It is basically a hierarchical model comprising many layers and subsystems. 


BUS-TYPE 
TL-SEND/- 
TL-RECEIVE 
TL-RESIS 
TL-REAC 
TL-SU SCEP 
TL-COND 
N-TL 

N-BUS 
N-SH 
SH-BUS|} 
SH-RESIS 


~ Input Dynamic Simulator 
Data Power Flow Data Results 


Figure D.1: Main user-interface of the DSI-1 in the MATLAB-Simulink environment. 
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CLOTS IOICIOII© 


NET-POWER 


5) 
MAX ( |D-P-Q) 


NEWTON-RAPHSON-ALGORITHM 


Figure D.2: Schema of the Dynamic Simulator block in Fig. D.1. 
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Table D.1: Symbols and descriptions of DSI-1 in Figs. D.1 and D.2 


Symbol Description 
TL-SEND Transmission line (send bus) 
TL-RECEIVE Transmission line (receive bus) 
TL-RESIS Transmission line (resistance) 
TL-REAC Transmission line (reactance) 
TL-SUSCEP Transmission line (susceptance) 
TL-COND Transmission line (conductance) 
SH-BUS Shunt elements (number of buses) 
SH-RESIS Shunt elements (resistance) 
SH-REAC Shunt elements (reactance) 

N-TL Number of transmission lines 
N-BUS Number of buses 

N-SH Number of shunt elements 

TOL Tolerance of the calculation 
BUS-TYPE Bus type 

IT Number of iterations 

COND-NO Condition number of the Jacobian-Matrix 
D-P-Q Active-reactive power mismatch 
V-M-0 Initial voltage amplitude 

V-A-0 Initial voltage angle 

V-M Voltage amplitude 

V-A Voltage angle 

P-LOAD Load active power 

Q-LOAD Load reactive power 

P-GEN Generation active power 

Q-GEN Generation reactive power 

Y-RE Real part of the admittance matrix 
Y-IM Imaginary part of the admittance matrix 
P-NET Scheduled active power 

Q-NET Scheduled reactive power 
P-LOOS-TOT Total active power losses 
Q-LOSS-TOT Total reactive power losses 
P-SLACK Slack active power 

Q-SLACK Slack reactive power 
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Appendix E: Software implementation of DSI-2 


Here is the implementation of the DSI-2 which is a main TR control 
system. This simulator is implemented in the MATLAB-Simulink 


environment with a user-interface, as shown in Fig. E.1. 


‘Ss Error 


VSt-ref a — 
CONTROLER & ACTUATOR TRANSFORMER 


RESULTS PLOT 


VS1-m V51 4 1S1-- — 


LINE-DROP COMPENSATOR 


Figure E.1: Main user-interface of the DSI-2 in the MATLAB-Simulink environment. 


Table E.1: Symbols and descriptions of DSI-2 in Fig. E.1 


Symbol Description 
PT Transformer active power load 
QT Transformer reactive power load 
VSO Primary transformer voltage 
VS1-ref Reference voltage 
VS1 Secondary transformer voltage 
IS1 Secondary transformer current 
VSi-m Measuring voltage after line-drop compensator voltage 
Delta-V Voltage error 
Tap Tap position 
a Transformer tap-ratio 
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Flexible Optimal Operations of Energy Supply 
Networks 


This book presents a systematic study consisting of modeling, 
simulation, and) optimization of dynamic operations of energy supply 
networks with distributed generation (DG) and battery storage systems 
(BSSs). Based on complex power flow models, different optimization 
problems are mathematically formulated and solved, In addition, novel 
mathematical models and a new combined problem formulation for 
active-reactive optimal power flow (A-R-OPF) in passive distribution 
networks (PDNs) (without DG units and BSSs) and! active distribution 
networks (ADNs) (with DG units and BSSs) are studied: Typically, 
distribution networks (DNs) consist of two different networks in terms 
of voltage levels, namely, low-voltage and medium-voltage DNs. For 
this reason, investigations are carried’ out separately on both networks. 
Modeling procedures for PDNs, ADNs, and energy prices are presented. 
These procedures serve as the basis for this work. 
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